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ABSTBAqT 

The application of latent trait theory, to classroott 
tests necessitate^ the use of siall .sample sizes for parameter 
estjlmation. Computer genet^ted data were used to ds^ess ^he accuracy, 
of estimation of the slope and location parameters in the tiro 
parameter logistic sodel irit^h fixed abilities and varylng^^smiaXl ' 
sa-mflie sizes. JcVe maximum likelihood procedure for estimating the 
parameters ira's compared to a method .in irhich the observed relative 
.fpeguenciesT irere smooth^ed using an isotonic regression method prior 
to applying, the maximum likelihood procedure. The isotonic method. ira&- 
considered 4)romising becausey^he^ smoothed ^el^^tive fregu^ncies" yield 
sore accurate ^estimates of €he probability of ccrrrectiy ansvering a 
test item'givena particular level of the ability than do the 
observed ^relativte .frequencies. The results were presented in terms of' 
variance and mean squared error of estimating the parameters. The 
results indicated that the' isotonic procedure provided, more accurate 
estimates of. the locsitiTO parameter vhere|is the maximum likelihood 
procedure pi/ovdded mor^ accurate estimates of the jfelope parameter. 
Since the isotonic method did provide for more accufate , -estimation of 
the location parameter # it^vas concluded that 'the isotonic method 
warranty further attention* The implications of the results for use 
irith classrobn tests irere also disoussedt (Author/CTM) 
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ACOIRACy OF. ESTIMATIS6 TS^HPARAMSTER XOGISTIC 



LAimr ISAfT PARAMEmS ANDv"lMPLICATIONS FOR CLASSROOM TESTS 

Michael J.' Kolen Douglas *R.' Whitney 
UnlTerglty o f _Ipwa 

Abati^act 



^ The application' of l^entVtifalt theory to ^dlaaaroom testa aegje^sftateb.^ 
the/use of sma3,l sataiple sizes for paraineter* estiinatl^r" - (^^^ 
Jata were used to assess the accuracy of esti^'^3U>ii 'oZ the .slope and^ location 

^ ; ■ - -J S 

parameter's la the two par^eter loglst^c^iodel with fixed abilities andr^vaVy- 
ing small .sample s^zes. The maxliom likelihood (H^ .procedure for eatlmatlug 
the parameters was compared tro a method In which tH^ observed relative fre* 
'quencles were smoothed using aa Isotonic regression method prior to applying . 
the HL procedure. Ihe Isotonic method was considered promising because the , 

smoothed relatly^ frequencies yield jEnore accurate estimates of the probability 

/. > ^ ^/ . • ■ ^ ■ , 

of correctly answering a test item given a parjticulaif level of ability then 

4o the observed Tj^latlve frequencies. 

' ^* ' ■ ' , * 

^> ' % ' 

The results. 'Were presented in terms of variance and mean^ squared error, of 

estimating the parametersv ' The results . in<|lcated that the^lS()tonlc procedure 

provided more^ accurate estimates of the location parameter where&s tt^e ML pro- ' 

cedure provided more- accurate ^estimates of the alope par'ameter* Slncfe the""' 

. , \ ' ' ^ , * ' \ 

isQ^onlc iqethod did provide for mor^ accurate estioiation of the location param- 
eterf- it was concluded that the isotonic method ^arranta further attention., 
Ste. in^licatlona of the r^eaults for u^e with' classroom testes were also discussed. 



^ ACCURACY OF ESTJMATING TWO PAKAMETER LOGISTIC, 

. ' • '*.'■ ^ ^ 

LATENT TRAIT PARAMETERS -AJfl) IMPLICATIONS FOR CLASSROOM TESTS 



Latent trait theory is an alternative to clas^sical -theory for evalu^ 
at^g classroom tesj:3* Receilt application^ and suggested applicatiobs- 
"of-iatent tfatt-jtheory (e'.g. Baker, 1977; Hambleton and Cook, 1977; . 

Lord, 1^77; Marco^ 1977; and Rene2f and^ Bashaw, 1977) Uave niQSt'often been 

" \' 

in the context of standardized tedt's or in other 'situations in whic'h a 

^ • ^ ^ - \ ^ ' ' r ^ . ' ^ : 

large nuinbef of 'individuals are'abl^e-to ta^e the te^t items* Since class- 
room tests are ^uite often administered to a 'small numSer of individuals, ' 
the problem of 'possible low accuracy^^in estimation of the latent tr^it 

parameters due to small sample sizes must &e considered before latent' 

^ ■ , > . 

trait theory can be usefql for classroom tests* 

* . , ' V . ' 

This^study was an attempt to assess, ^using computer generate data, 

* *• ' . ■ 
the accuracy of estimation in the two parameter logistic model with yar^^ 

ing small sample sizes and to ex^^e one methad* offering ptomise of^ 

. , ^-t f ' . ' \ ' ' ' 

improvingPthe accuracy* Specifically,' £he frequently used maximum like-^ 

' ■ i ■ ' ^ x^* < ^- • 

llhoaid (MLX procedure for estimatijiig the /^aterft trait parameters was 

compared to a method In which the observed relative ft-equencles were 

smoothed using an isotohic regression metho4 l^rior to' applying the ML 



procedure;^ In the i^dtoriic method, the smo^hed'Telative frequencies (instead 
of the observed relative Trequencies) were used^^ the dattf'in the usuajl 
maximum likelihood. procedure^ ^ The Isotonic m£rthod was-^ considered promising^ 

because it results in^ the smo^tb^d- relative freqUenciesVwhlbK are entered 

' ' ' - . " , " ' ' ' "\ / . ^ I ' ' 

into the maximum likelihood prqc^dure' to^^be monotonlc non-decreasitig - 



along the ability j&cale* 



■■■■ r>- ■: 



It was believed that the Is'otonlc method ittlght le^d to tiiore accurate 
estimation of the latent trait parameters because smoothed relative fre- 
quencies exhibit a monotonlc nondeGreasing relationship that Is consistent 
T^ltji^he assmaptlons of the ML procedure* !rhey might, therefore, yield moye 
accurate estimates of the actual prdbablli^es than do the observed* re^lative 
frequencies (Ayer^ Brunk, Earing/ Reed, and Silverman, 19^5)1 l^e ysd of the 
smoothed relative fr&quencl^s also ensures^that a non-negative slope estlidate 
'will result when the. logistic paifameters are estimated « ' : - * * 

Comparisons of accuracy in estimatlbn were based on meaa squared / 
error, variance, and bias in estij^ting the parameters* In addition, the 
obtained average variances were compared to an approximation to the Rao-Cramer 
lower bound £or the variance in estimating &ach of the two parameters* 

. . METHOD 'AND J)ATA GENERATION 



This stu^y uded the two parameter logistic tnodel (Blrnbaunit 196d)~:' 
exp U:? a (0 T b )} 



g 



1 + exp U.7a (0 - b )] 



(g 



' > ^ 2 , * * * , It) 



(1) 



In this model, P (0) Is the probability that; an examinjee with ability 0 
answers item g correctly^ a and b are parameters asso^ated with item 
g (g * 1^ 2, k), and,k is the numbei: of items on ihe test* The 

parameter a (often tefecred^to as an index of item dlsc|rlmltl^tloil) is 

I g ' * * * 

a slope parameter and b i>ften referred to' as an index ol^ item, dlfflctilty 



representing the ability level at which 50% of the individuals Iwould b^ 
expected to correctly answer the item Baker# .1977) is a loaatiorf f 



parameter. The constant of 1.7*'waj^ chosien so that the probabilities and 
paratneters s^oulH'cprrespond closely to those of a noma! oglve^ 

Marginal dlstglbiit'lon of ability . The marginal distribution ^ * 
Ability was assumed 'to be a logistic distribution with the distribution- 



function : 



exp (1*70) 



jr^mt^^^^pl^j^ was divided into 15 equal probability ,intervals 

^that were 'fixed for all replications* The median of the interval was 
used as th^ ^ability level of all scores in that interval^in subsequent 

' calculations* ' * - 

\ ' \ . * ' ' * 

Item c haracteristic curves: The actual item characteristic curves 

were conatruc^d'using the two parameter logistic models with a^. taking 

the values 0^5, 1.0» and 1*5 and b "taking the values 0*0» 0.5» and 1*5. 

Sample size (S) was allowed to £ake on the values 45^ 60,*90, 150, and 

' ' * , d t ; . \ 

.240^^ These values were chosen td cover the range of possible item para- 
meters and sample sizes that are likely to occur with classroom tests. ' 
/ , Random, generation of the relative frequencies ^ After the logistic 

'paramet^ers and^ sample size were fixed* the relative frequency of individuals 

' " ' w» / ^ * ' 

^ j:o];r€fctly answering the ,item within each of the fifteen intervals on the' 

" kbHity scale ^was ^ranct^mly generated using-the IMSL (1977) Fortran IV ' , 

subroutine/ GCBLJt^f ^ The generated binomial variates were based cfn*the 

actual P (0.0 value^£jrom the two parameter logistic ^model (i*e*» the 

probabilit}^ of cor^ctly Answering the* item for individuals ,^t the median 

o*f the iriterval)>ana^:5th^? satople size in the interval* The randomly 

geneifaEed 'yaiue.fcprrasp^c^nded' to the relative freqixency of Individual^ * 

with ability*' 0 that ^perfectly answered the hypothetical test item. 



. Maximum likelihood eyimatlon methods . Maximum llkellhqM procedures 

were used to estimate th^ a and b parameters* The llkellhoad function 

for the two parameter loglsjH.<v model for^an individual item is prcfportiodal 

to e^ (Finney, 1971) where, 

m m f ^ 

L r 'In [P-^ (9.)] + I (n r . ) -In [1 - P (0 )], 

1 represents onb of the m = 15 -intervals along the a}>ili"ty scale, 6^ J^* 

' % ' s ■ a 

the median Ability for the .th interval, n. Is the' numbelt'^of individuals 

in interval i, r.^is the number answering the tteifi correctly, and P 

" It . ^ , _ * * ^ 8 

Oj) is as given in equation/(l) . The ma^tlmum likelihood solution 

ijivolves finding the a_ and ^ values which maximize L. maximum 

% %j — . \ 

likelih<^d equations $re obtained^by finding the partial derivatives of 
L with respect to a and Jb and setting them equal to "zero. ' Direct ^olu-, 
tions are not possible, so the values of a and b which maxlmize-L ^rc 
found by iteratively solving these two maxiimim likelihood equations; 



2 0.[P^(0.) --^1 = 0 



(2) 



2 [P (Q.) - T^l 



1=1 



:8' 1 



The ntaxisnim likelihood 'solution was acftleved by substituting stantjtiig*' 



values for and b Into equation (1), the equation for (0^),^^^ajw,, using 
an Iterative scheme to find the«it and 6 (the edtlmated 'phram^et) <calues . 
which maximize L, * ^ * ' ' 



parameters a :and b were estimateud using the IMSL^ (1977) Fortran IV 
. ' S g ^ . ■ / 

subroutine ZSYSTH which uses Brown^s algorithm*^ The convergence criterion 

■ ' ' ' ' / ' 

vas satisfied lf» on two succefss^ve iterations* both estimates' agreed to 

^significant digits* Ip order to minimize the number of iterations and 

'lack of convergence problei^s the actual/yalues of a and b. were i^ed as 

g p S , 

the starting values for the Iterative process. Fortran IV double precl- 

• '' ^ . • ' 

slon arithmetic was used* for all of the computations* , 

\ ■ • ^ 

Parabe'ter esfTTnation itethods * Twp procedures were used to estimate 
the param&f^s* In the first (regu),ar)» the masc^um .likelihood solution 
was applied dlrfe<;tly to the observed relative frequencies. In^the second " 
p(^ocedure (isotonic) the ohBetwj^ relative frequencies were dnltially. 
smo</thed using an isotonic* re^&ess^n method* In this procedure* the ran- 
domly generated ^relative fi;equ^cles oi. individual^ correct^ answers were 
checked to determine if ttfey were monotonlc nondecreaslng along the ablllt>^, 
scale* Whenever a pair of relative frequencies did not meet this condi- 
tion (a reversal)/ the pair was averaged ^ The procedure continued until 
no reversals remained * The maximum likelihood solution wa^ then applied to 
the siip'othed relative frequencies. 

Data generation * The relative, frequencies were randomly generated 100 
times for each of the 45 possible combioatiops of a » b » and N* The 
parameters (a and b ) were estimated *for each of the replications using ^ 

both of the procedi^res* The' individual ability parameters were fixed and , 

* 

assumed known for all replications* 

Assessing the accuracy of estimation * The accuracy of estimation was 
judged by comparing^the mean squared errors (mean square'd devia^tions about 
the actual parameter values) and variances (mean squared deviations abotit ■ 
the average €^timated par^eter values) in estimating a and b for each 



combination of the logistic parameters^ sample slze» and estimation ^v^.^ 
method* Meati values of the parameter' estimates were also examined- in or<ler 
to ^^ess the amount of bias resulting from each of the esti^i^ion 
methoils; The Xpglstlc parameters on some of the replications were not ^ 
estimable by the maximum likelihood method*^ For comparative purposes^ ' / 
only those xeplications for which both the regular and isojtonic methods ' 

^proved estist^ble were included in the computations of the summary stSatis**- 
tics, ' ' . , ' 

Approximating variances* the formulas for^which are derived in the 
Appendix* were also calculated. 'The approximating variance for the slope 

' parameter is given by - ' ^ ^ r 

2 ° 2-1 
Approx/ Var. (a ) = [1.7^1 w (0 - b^) ] 

& r 1*1 ^ ^ O ' _ 

and the approximating^ variance for the, locat;Lon parameter is given by 

Approx, Var, (b .) = [ (1*3^ a Z wj 
* - — * i*>l ■ 

where ^ 

These approximating variances hold vhen the P (6,) values follow 'the two 
parameter -logistic distribution* the samples are rtodom at. each ability 
level wi^ ability levels being fixed* and the number of individuals at^ - ' 
^ch abiXtcy" level is large (Berkson* 3953), Because the first two of these 



- .1 



3^gn 



conditions vere met in the design of the present study*^ the results ma/ 
be used to assess how large a sample stit^ is needed for the approMmating. 
variances to provide close approKlmations to the observed variances < 



> Results ■ 

' w ' ' ' * 

•A solution was attained for 99.36% pf th^ solutions atteiiq>ted using 

-the regular method and , for 99.93% of'the ■^oO^tio.ns attempted using^the . 

isotonic method. The replications in which. solutions were not found 

tended to be those with sznaller sample sizes and larger parameter/values*\ 

The only combination of parameters and sample size In which fewer than 

90^ of the solutions attempted were attained was with N=*45, a =1*5, and 

^ 8 

bg=1.5.. For the regular method, 81% were s^ved and for' the isotonic r 
method, 83% w^e solved in this case. Only, t^iose items for which both 
the regular and isotonic metjjods produced solutions were included'in sub- 



sequent analyses. For these items, the ayerage niim^er of iterations was 



82 fpr the regular method and 5.?2 for the isotonic method. 

Mean estimates of the parameters. * * The mean values of the slope 
pg^"i^^]p^'f>^r fa ) &cl^^maJKAla_ar:A---prpflpn^A/^ "in-Table^i and th^ mean v^ues of - ' 
the location parameter (b ) evstimates are presented in Table 2. Fox all 
combination^ of mett^od, sample size, ^and logistic parameters except^ne 
' (regular ^^thod, N=15(T, a ==0*5, and b *1*5), the mean estimate of the " 
'slope parameter was greater tlian the respective actual slope parameter* 
In addition, the isotonic method produced larger estintates of the slope 
parameter than did th^^regul'ar method in all cases* Ths bias evident, 
by examination of Table 1 is in Che same direction as the bias reported by 
Berkson (1955)* The^e was a weak tendency for the bias to lessen as 

* 

sample size increased* ' ' " ^ ' ^ 

The means, of the estimated location parameter values .do not appear to 
be consistently biased in either direction (Table 2).. However, the means for 
regular method vere consistently greater in absolute value than were the . 
means for the isotonic method* 



Insert Tables 1 ^nd 2 about here* 
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Accuracy of Estimation 



r 



The variance and approximfl^g variance for the slope and locEftlon par'a- 

' . , - • ' . ■ - • ; I ' 

'meter esti^mates are presented In Tables 3. and 4. The mea^ squared ^errors 
In the slope and location parameter estjjnates -a^e presented In Tables . - ' 

5 and 6i ^ ' ' , ■ . / " I ' ' 



Insert Tables- 3 » 4» S» and^ 6 about here 
-w— 1 ; — ^ — 



^' Slope garaiPeter * The regular method provided more accurate^'esTtimates 
of the glope parameter than dl<^the Isotonic inethod (Tables 3 and '5) for 
most combliiatlons of tjie paramet^s and sample &lzea« The relatlVe 

differences lb accuracy for the two methods were greater In mean squared 

• / 
error than In' variance prlraazlly because the isotonic estimates of the 

' ■ s " . . ■ . ■ 

slope parameter were more.blased than were the regular estimates* The 

relative difference Ih ^ccurady for the two^metW^ decreased as sample 

size Increased* (This Is as expected^ slnce» as the sample size Increases* 

-^ne observe'd relative frequencies less frequently require smoothing*) 



/ 



location Parameter , The isotonic method ^provided more accurate 

estimates o£ the location parameter than did the regular method (Tables 

* 

4 and 6) for all combinations of the parameters and sample 'sizes* The 
A . * o . ■ 

relative differences between the mean squared ^errors and variances were * 
in general, very small* 

' \ Approximating Variances * The variances for slope parameter estimates 
were in all caaes larger than the asymptotic variances (Table fS* The 
approximating variances substantially underestimatl&d the obtained variances 
for smaller ^samples (as\:exprecte^) and larger^ parameter values * The 
variances (and mean squareci. errora) for the isotonic method wer^ very 
sUmiXar to approximating ifariances for Estimating the location parameter. 
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The * Isotonic me Alod* tended to r^^uZt. In. variances that were actualJby ^ 
lowet than the approxlioatlng va^r^ancee for Small samples and lower values of 
the location parameter/ The«^ variances' (and'pean s^iuared errors) for the 
regular method were similar to t^e approximating variances for the larger ' ' 
sample sizes* ■■ * - ^ . . \ . / ^ 

Di^coission - , 

• . ■ ■ ■ ■ \ ; - ' . . 

The ayejage slope parameter- estimates vere larg;^r 'for the isotonic *^ 
'^'method than for the regular -method « 'A possible e^latiati&n for this 

result may be' attempted using an analogy to linear least squares regr^ssiofi 

■ */ ^ 
.analysis* The latent trait maximum likelihood estimation problem may Be 

solved using weighted linear regression analysis, with the' predictor being 

> 

ability and the criterion being 1,7a (O^-S ), In this solution, it is 

^ ^ ' ' . ■ ' ^ ' ^- , / ' ' 

necessary to iter^tively splve for the parameter^ and fot the weights which 

'axe/a^^fimcfrLon of the parameters (for example, see Befkson, 1933)* The 

* - ■ t 

isotonic tnethod initially results in. the predictor and criterion being in 

the same nondecreasing monotonic ordering due to averaging the criterion values 

involved in reversals. If the isotonic'method was used prior to linear 

l^st squares regression, th& t:Qvariance between predictor and criteri(fn 

(and, therefore, the slop^) would increase in compari^a^to that which w 

would be otM:ained \f tixe isotonic m^hod had not bedlHsed, (The variance 

of the predictor remains udctianged^) If this analogy^ -^"^d reasoning holds, 

it would e^cplain why the slope estimate is*^eater under the isotonic methods 

l^e smaller location parameter values foun(} with the isotonic method 

a^e more difficult to e^cplai^ rsiay be that the slo|^e*and location ^ ^ 

parameters are related in the iterative scheme so that an increase in the 

slope result^ in a dgc^re^^^ in/the locatj.on parameter. 



- The regular method led to mo]re acctirat^ e^tima^^loi^ of the slope 

parameter where a's the Isotonic mfetljg^yi^lded more accurate estimation 

^ t ' ' ' • , ^ 

of ^be location par^eter. The Isotonic method was no^ more acfur^te 

for ihoth .paramet^9 and does not provide a reduction in, comput^l^lonal 

labor (no iFeduct'l-on was noted In the nttmh^r iterations required) \ 

' \.- . ■ .. ■ • ; 

It probably, ^"therefore, should not .he used .exclusively in practice, ^ ^ 

*A^aieme that would Te heneflcial would he to fit the data with both 
methods inil to use the regular tiuethod to ^estimate the ^siope parameter 
and the isotonic method to estimate the location parameter.* However, 

%he gain in accuracjr over using jtist the regular method might hot 

justify the ^UJdltional' computational labor. Another scheme would be ^ 

[ * 

tot estimate the parameters using the isotonic meth^. Ignoring the 

.isotonic estimate of the slope, the isotonTc ^estimate qf "th?e locatS^ 
parameter wc^ld then be used as a fixefl value for estimating the slope 
using the observed relative ftequencies as data. 

* 

''The results also indicated that the approximating variances approach 

closely the obtained variances for the larger -Sample sizes and at some 

combinations of the parameters for the smaller ^Ample slzes^* Thus, the 
• ■ "■^.'^^^v - 

approxlmatlng variartces can be used to indicate the degree of accuracy to 

^ \ \ ■ . / 

be expected when using two parameter latent trait ^theory with fairly 

large' sample sizes^ The ,use of the aj>proxlmating variances assumes know- 

ledge ^£ stuiafent ability, that the logistic model actually holdsj and 

* ' ' • .... ' ■ 

/that there is no^bia^ in estimation* These are strong assumptions to 

ma(ce in applied settings '(and in fact, the estimates of the slope will 

be biased)'^3o tha£ the approKimatln^ variances 'should be seen ^as a -lowef 

bdund which would probably not be achieve^ — eyen by the regular method, 

^ ■ ■ ■ ^ ^ ' 



In rel'Htlon to classroom testings unless modera^tje sainple sizes are 
used (say 100 or more students) the latent trait parameter estimates will 
'b^. fairly Inaccurate - (with; the degree Qf accuracy also depending on^the « 
actual parameter valuers) However » the problent of Inaccuracy In estimating 
the classical Indexes also'^xlsts (B^keri^ 1965),* ^In comparing the results 

of the present ^tudy to tho$e of Baker' (1965)» It appears that the estimates 

. ■ * •■^ 

of the latent trait parameters are probably no' less accurate than those 

for the classlc'al Indexes for small samples* Thus» If It is desired to 

present analyses to Instructors*' no clear preference emerges for 

another comparison of the relative accuracy of astimatlon In the classical 

* I. 
and latent trait theories* , ^e choice should be made on the amount and 

^ ^ \ ' ^ " ^ ' . ■ 

quality of Information each. theory provides to th^ Instructors* . 

' \ ■ " ' - ' ^ . " 
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^ Table 5* Mean squared 



error for the esjtimated slope (a ) parameter. 
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Table 6* Mean aqtiared error for the 'estimated location (b ) parameter' 
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Appendix: 



A metliod which provides an approximation to' the Rao-Cramer lower 

hound for the variances 'of & and ^ (the Item subscript; g» wlll.be 

8 8c.* " 



omitted) were used* 



Berkson's (1953) description of a proof for a s&ghtly different 
{farameterlzatlon of the logistic model was. followed* Plxst^ some results 



will be stated regaling the^ logistic distribution which will be used In 



< 



the derivation* ' As stated in the text^ the distribution function for the { 
logistic distribution used here is: * ' ' ^ 

P(e^) = exp [1,7a (Gj-^b)] /. {1 + exp U,7a (0^ - b)]}. ^ (1) 
If both numerator and denoifiloator are divided by exp [l*7a (Q. - b)'J then, 



P(0^) = (1 + exp [-1.7a (0^ - b)]} 



-1 



Also, 



(3) 



1 - P(0^) = 1 - exp [1.7a (0^ - b)) /^{l + exp [1.7a (0^ - b)]} 

1 + exp [1.7a O' - b)) - exp [1.7a (0 ' —b)) 
= . ± : i . 

, {1 + exp [X.7a (0^ - b)]} 

= tl + exp tl.j> (e -■b)jr^. 
ripall^, dividing (1) by (*3) we have, 

P(0j_) / U - P(e^)] = exp' {1,7a (9^ - b)]- (4) 
To find the approximating lower bound to the variance, it is neces- 
sary to examine the likelihood function. First, thia likelihood function 
'will be placed in a more convenient form, * . 
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From, the text, 



L = Z r,-la'pf0.) +1: (a, - r , ) *ln [Iv- p'(G,)] 
1=1 ^ ^ 1=1 ^ ^ ^ 

m ''ill m ' 

= I r.-ln P(0,).+ I n. In [1 -P (0,)] - I r.'ln tl - P (0,)]. 
1=1 ^ ^ 1=1 ^ ^ 1=1 ^ ^ 

Combining the first and third terms, 



-1 



L- r r.'ln {PCe.) / [1 - P + I n. In [1 --p -"(e,)] ' 

. 1=1 ^ _ - ^ ,1=1 ^ ^ 

Sy substituting (3) and (4) J.nto the equation, ' ■ . 

.L = I r.-lh {exp [1.7a (0 b)]} + I n lii {l + exp [i.7a(0,rb)]} 
1=1 ^ ^ • 1=1 ^ ' , ^ . 

=. 1.7 Z r, a(0, - b) - . S n. In {l + exp tl'7a (0, 7*b)]}. (5) 
1=1 ^ / 1=1 ^ ^ 

To obtain the. lower bound for the variance of 3, It Is necessary to 

find the first and second partial derivatives' with respect to a. Taking 

*' 

the flrst'^derlvatlve we gfet, • . ~ 

37 ' m . _ 

= 1.7 I r. (0, - b) - 1.7 £ n, (0, - b) {l + exp { 1.7a (&.-b)J} 
3a 11 j_ .1 1 1 



• exp [1.7a (0. - b)f. 
By using result <2) and rearranging terms, ' • ^ , 

1^-1.7 Z_ r^ (0^ - b) - 1.7 I .n, (0, - b) {1 + exp [-1.7a<0,-b),]}"^. 



l..^"l 



1-1 * " . 1=1 

Taldng the second partial derivative with respect to a. 



0-1.7 Z n, (B. - b) (-1) {1 + exp {-1.7a (0..- b)']}"^ 
3a 1=1 ^ ^ 



• exp [-1^7a (0^ - b)]; • (-1.7).(0^ - b)- 



*- -(1.7) t n^ {0^ - b) • . {1 + exp |-1.7a (0^ - b)]} 



-2 



1-1 



23 *' exp {-1.7a^(0^ - b)]. ' 



- - 21 - 



Note that* 



{1 + exp [-1.7a (0^ - b)]}"^-- exp [-1.7a {0 - b) }, 
- {1 f exp [-1.7a (0^ - b)]}"^ • {1 + exp [-1.7a (9 " - b)M exp [i-.7a '{Q, 



,-1 



-1 



• {1 + exp [-1.7a (0^ - b)]}_ • {1 + exp [1.7a (0^ -b)]} 
^ow, using (2) and (3) the above expression reduces to^.^ 

• pxep [1 - p (0^1. ■ . . 

' • ^ . - ' 2 ' 2 

Therefore, substituting back Into ^he equation for 3 L/3a we have. 



(6)" 



3a i^l 



Define, -w^ - P(e^) [1 - P (0^)]. ^ 
Then, , . " 



9? 



m 



-i. 1(1.7)2 i (0^-b )2 
3a 1=1 

and the approximate lower bound [which 
estimating a ls» ^ 



2 ^ ^ 2 ' 

(1.7)^ Z w, (9, r b„)^ 

1-1 * 



], for the variance In. 



(7) 



To obtain the approxlmat:e low^^b'ound 'for bj It Is necessary to find the 
first and second partial derivatives of L with respect to b. Using the form 
of L In equation (5), 

ftT mm' 

^ f -1.7a 'I r , - I n, {1 + exp [1.7a ($- - b)]} exp {1.7a (0. - b)] • (-1.7)a 
. 3b 1=1 ^ 1=1 ^ ' ^ 



By rearranging terms and using results (1) and (2), 



- -1.7a Z r^ + 1..7a Z n^ {1 + exp^[-1.79 (0^^ - b)]} 



-1 



1-1 



1-1 
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0< 

Finding the second partial with respect to b» 



^ - 0 + 1,7a S n. (-1) {l + Bxp [-i,7a (G, - b)ir^ 
3b^ i=l ^ ^ ^ ^ 



Rearranging terms » 



1^ 



. ■= -(1.7^)^ S n. P(0.) [1 - P (0,)]. , * ; 

By taUng,|w^ = P(e^) U - P(0^)], 

3b 1=1 ^ ■ 

Then the approximate lower bound to the variance In estimating b Is, 

1 • 

2 

1=1 ^ . ' 



exp [-1.7a (0 
• (1.7a) 
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